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^ , Abstract 

^vq \ We apply quantum-mechanical sum rules to pairs of one-dimensional systems defined by potential 



energy functions related by parity. Specifically, we consider symmetric potentials, V{x) = V{—x), 
and their parity-restricted partners, ones with V{x) but defined only on the positive half-line. We 
extend recent discussions of sum rules for the quantum bouncer by considering the parity- extended 
version of this problem, defined by the symmetric linear potential, V{z) = F\z\, and find new 
classes of constraints on the zeros of the Airy function, Ai{C,), and its derivative, Ai'{C,). We also 
consider the parity-restricted version of the harmonic oscillator and find completely new classes of 
mathematical relations, unrelated to those of the ordinary oscillator problem. These two soluble 
quantum-mechanical systems defined by power-law potentials provide examples of how the form of 
the potential (both parity and continuity properties) affect the convergence of quantum-mechanical 
sum rules. We also discuss semi-classical predictions for expectation values and the Stark effect 
for these systems. 
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I. INTRODUCTION 

Sum rules, constraints on energy-difference weigfited combinations of on- and off-diagonal 
matrix elements, have been an important theoretical construct in quantum mechanics since 
its earliest days. For example, the Thomas-Reiche-Kuhn (TRK) [l| sum rule was one of 
the earliest quantitative checks on the oscillator strengths of atomic transitions and an 
important confirmation of the applicability of quantum mechanics to atomic systems. Their 
use continues today, for example, with QCD sum rules used to probe the masses of both 
light and heavy quarks in hadronic systems. For a list of recent applications of sum rules, 
see Ref. |2|. 
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The study of sum rules and how such constraints are realized [2j can provide useful 
examples of a variety of mathematical techniques utilized for their confirmation. In addition, 
sum rules can be used to generate new mathematical constraints [Sj on quantities involving 
infinite sums. For example, in exactly that context, two of us have recently shown how it is 
possible to systematically derive new constraints on the zeros, (— Cn), of the Airy function, 
Ai{(), of the form Sp{n) = J2k^niCk — (n)~^ for natural values of p > 2 using quantum- 
mechanical sum rules applied to the so-called 'quantum bouncer' problem |J]. Defined by 
the potential 

{Fz for z > 
(1) 
oo for 2 < 

this quantum-mechanical problem has been applied to recent novel experiments, ranging 
from the bound states of neutrons in the Earth's gravitational field J5| to optical analogs of 
gravitational wave packets |6j] which demonstrate the presence of predicted 7| wave packet 
revivals. 

In our earlier study, a number of the most well-known sum rules did not converge due 



to the lack of continuity of the bouncer potential at the origin. In this note, we extend 
the results of Ref. [4[ to discuss mathematical constraints which arise from the study of 
quantum-mechanical sum rules in a closely related system, namely the symmetric linear 
potential, defined by 

V{z) = F\z\. (2) 

This is an example of a parity- extended potential, namely taking one that is defined only 
for ^ > 0, and extending it in a symmetric way over the real axis. In this case, we find new 



relationships involving both the zeros of Ai{() and its derivative Ai'{(), namely —(n and 
—rjn, respectively. In Sec. Illll we also explore the convergence properties of sum rules in this 
system, compared to the quantum bouncer, because of the improved mathematical behavior 
of the potential energy function in this case. 

For comparison, we also consider in Sec. |V] a parity-restricted version of the most well- 
studied system in quantum mechanics, namely the harmonic oscillator. Specifically, we 
examine the structure of sum rules in the 'half '-SHO potential, defined by 

{muj^x'^/2 for a; > 
(3) 
for X < 

and find dramatically different behavior for many sum rules. 

We are motivated to study the relation between such parity-related potentials since the 
solution space of the parity-restricted version of a symmetric potential, V{x) = V{—x), 
consists solely of the odd energy eigenvalues and eigenstates, with only a trivial change in 
normalization. Thus, half of the wavefunctions and energies will be closely related in the 
two systems. Since sum rules involve intricate relationships between energy differences and 
matrix elements derived from the wavefunctions, it is instructive to see how identical sum 
rules are realized in two parity-related systems. We find that many matrix elements can be 
very different in character between two such related systems, yielding sum rule identities 
which are realized in novel ways. 

For these two systems, we then make contact with semi-classical (WKB-like) probability 
distributions (in Sec. IVI|) and examine the similarities with the new exact quantum results 
presented here in the large n limit. We also note that a recent examination of the Stark 
effect in the symmetric linear potential [8|| has found interesting constraints on second- 
order perturbation theory sums which are very similar to the energy-weighted sum rules we 
will discuss below, and we extend that analysis to the 'half'-SHO system introduced here, 
allowing for insight into the structure of the Stark effect in that system. We begin by briefly 
reviewing the results of Ref. 8| in the next section, mostly to establish notation. 

II. SOLUTIONS FOR THE SYMMETRIC LINEAR POTENTIAL 

We begin by reviewing the properties of the solutions for the symmetric linear potential 
in Eqn. ([2]), using a slight modification of the notation in Ref. |8|]. We note that the only dif- 



ference in notation with Ref. [^ is in the normahzation of the wavefunctions, corresponding 
only to a physically irrelevant difference in phase. 

Using the fact that the linear potential admits Airy function solutions, with different 
forms for positive and negative values of position, we can construct piecewise-continuous 
solutions of the corresponding Schrodinger equation. Because of the symmetry properties 
of the potential, we can also classify the solutions by their parity, with the odd functions, 
which must vanish at the origin, related to the solutions of the quantum bouncer problem. 

n n 

For background on the solutions to the bouncer problem, see also Refs. [7[, [9[ 
Specifically, the odd solutions can be written in the form 



io|. 



^. 



(-)( 



Ni-^ 



(4) 



Ai{z/p - /3n) for < z 

-Ai{-z/p - /3n) for z < 

where p = (h'^/2mFy^^. The appropriate boundary condition at the origin is that the 
wavefunction must vanish there which implies that 



^(-)(0)=Az(-/3„) 



(5) 



where we identify /3„ = (n where the —(n are the zeros of Ai{(). The corresponding combi- 
nation for even parity solutions is written as 



i^i^\x) = N^-^^ 



Ai{z/p - I3n) 



for < z 
for z < Q 



(6) 



^ Ai{-zlp - /3„) 

which at the origin must satisfy [■?/'« (-2 = 0)]' = 0, implying that /3„ = r]n where the —7]n are 
the zeros of Ai'{C,). 



Handbook results ll|], 12| give approximations for these zeros, in the large n limit, as 



Svr 



(2n - 1/2) 



-1 2/3 



and 



37r 



((2n-l)-l/2) 



-1 2/3 



(7) 



where the labeling starts with n = 1 in both cases. We note that WKB quantization gives 
the same results J8|, valid for large quantum numbers. 

The normalizations required in Eqns. (jl]) and (|6]) are obtained by using integrals first 
derived by Gordon |l13| and Albright [ij] collected in the Appendix in Sec. |Al specifically 



Eqn. (IA7|) . and are given by 

Ar(-) = 



1 



2pAl'{-Cn) 



and 



iv(+) 



y/2prin Ai{-rin) ' 



(8) 



We note that the wavefunctions at the origin satisfy 

V'(-)(0) = and ^pM(0) = -±=. (9) 

The energy eigenvalues are then given directly in terms of the (n, Vn by 

Ei-^=£oCn and E^ = ^o^^n (10) 

where Sq = pF and the energy spectrum satisfies 

Ej+) < e[-^ < E^+^ < E^^-^ <■■■ . (11) 



Just as for the bouncer system pi| , we have a power-law potential of the form V(^k) (x) = 
Vo|x/a|^, and the quantum- mechanical virial theorem, (T) = |(V^), is satisfied for A; = 1 
which we can confirm by direct evaluation. For example, for odd states we have the matrix 

elements 

"+00 



ii^i^lVizMi^) = / \i;i'\z)\'F\z\dz 



2Fp 



2 



2p[Az'(-C„)]2 io 



y[At{y-Cn)Ydy 



= '^SoCn = Ie^^ , (12) 

using the integral in Eqn. ( JA8I1 . with a similar result for the expectation value of the potential 
energy for the even states. The expectation value of the kinetic energy operator, 

can be evaluated by either: (i) using the definition of the Airy differential equation, Ai"{(^ — 
P) = {( — l3)Ai{(), to rewrite the integral in terms of the one in Eqn. (lASp or (ii) using an 
integration by parts to move one derivative onto the first wavefunction, and then using the 
integral in Eqn. (lAlOp . 



The dipole matrix elements needed for many of the best-known sum rules have a different 
form than in the quantum bouncer. Using the symmetry of the potential, we have 

(^i-^kl^i-^) = (^i+^kl^i+^) = o, (14) 

while using Eqn. (IA12p . we find that 



As with all of the matrix element expressions and sum rule identities we present, we have 
confirmed these expressions numerically using Mathematica® . 

In contrast to the quantum bouncer problem [j] where all off-diagonal matrix elements 
gave terms proportional to inverse powers of (Cn — Cfc), using sum rules involving only dipole- 
matrix elements for the symmetric linear potential will give constraints on sums of powers 
of l/{i]k — Cn)- Matrix elements for even powers of z, for example 



¥r^V\i^^rt') 



(16) 



will connect even-even and odd-odd states. The corresponding odd-odd sum rules will 

n 

reproduce some of the same identities found in Ref . [^ , while the even-even sum rules will 
generate new identities, involving sums over inverse powers of {rjk — rjn)- 



III. NEW CONSTRAINTS FROM THE SYMMETRIC LINEAR POTENTIAL 



To illustrate the range of new constraints which the symmetric linear potential places 
on the (n and rjn, we first consider the most well-known sum rules involving dipole matrix 
elements. The famous Thomas-Reiche-Kuhn (TRK) sum rule pj] is given by 






n\z 



2m 



(17) 



while a completeness relation for the matrix elements of the momentum operator can be 
written in terms of dipole matrix matrix elements as 






ln\f\n) = n:- [E„ - {n\V{z)\n)] . (18) 



We note that since the states n and k that contribute to the sum rule will be of opposite 
parity, the k ^ n restriction is unnecessary and we can sum over all intermediate state k 
labels. 

Two other sum rules involving dipole matrix elements have been discussed by Bethe and 
Jackiw [15|, [l6|, namely 

Y,{E,-Er^f\{n\z\k)\^ = ^ln ^'^^"^ 

k 

and 



dz^ 



n 



(19) 



Y,{Eu-Er,f\{n\z 



h" 



m^ 



n 



dV{z) 
dz 



n 



(20) 



which are sometimes called the force times momentum and force-squared sum rules, respec- 
tively. We note that not all such sum rules necessarily lead to convergent expressions. In 
Ref. [^, Eqns. (TT9|) and (l20l) could not be applied in the quantum bouncer system due to 
the discontinuous nature of the potential in Eqn. ([T]). In the case of the symmetric lin- 
ear potential, the system is defined by a continuous potential and even the discontinuous 
derivative will lead to a convergent result in the corresponding sum rule in Eqn. fITI?]) . These 
results regarding convergence are consistent with the form of the dipole matrix elements in 
Eqn. ( !T5|) . where their values for fixed k or n clearly decrease more quickly with increasing 
n or k than their counterparts in the bouncer system, indicating better convergence. 

When we fix the state labeled n as being one of the odd solutions and sum over k values 
corresponding to the even states, the l/^yffk factors remain inside the summation, giving 
sums of the form 

rM^Y.^^^. (21) 

In contrast, when the even states are fixed, switching the roles of n and k in Eqn. (fT5|) . 
then the common l/y/rj^ factor can be removed from inside the summation and eventually 
moved to the right-hand side of the identity involved. In those cases, we find expressions of 
the form 

f^- W - E w.L^P = (^) E 77-^ - -UAn) (22) 



allfc ^"(^^ ~ ^'^^^ ^^"^ flit "^^^ ~ ^"^^ ^'^ 



where 



C^-W-EtT^- (23) 

allfc ^^f" ^"^ 

The related quantity where we sum over k values corresponding to odd states, namely 



allfc 

can be obtained by using the definition in Eqn. ( 12T|) to write 



w-Ei^- (24) 



Tp{n) = (n Tp{n) + Tp_i(n) . (25) 
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For completeness, we recall from Ref. (^ that for the quantum bouncer system constraints 
on sums were of the form 

^^-EtT-^Tv- (26) 



fc^n (^'^ ~ ^" 



We will focus on evaluation of the Tp{n) and Up{n), noting that it is straightforward to 
obtain the related Tp{n), Up{n) for comparison to Sp{n), using Eqns. fl22l) and fl25|) . 

While not formally a sum rule, we note that the completeness relation for dipole matrix 
elements can be written as 

Y,\{n\z\k)\^ + {n\z\n) = ^ |(n|2|A;)|2 = {n\z^\n) (27) 

ky^n all A; 

and that because of the parity of the potential the 'diagonal' term is in fact zero. While not 
a sum rule in the classic sense, it will also provide novel constraints on the Ai and Ai' zeros. 

For the sum rules in Eqns. (fTTl) - (!20|) and (1271) . we can now fix a state of definite parity, 
either ^n (z) or ^„ (z), and then sum over the states of opposite parity, namely the -ipl (z) 
or ipl. (z), respectively. In this way, we obtain two different constraints for each quantum- 
mechanical sum rule. 

For example, starting with the TRK sum rule in Eqn. (ITTl) we find the relations 

by using odd {^n ) and even (i/ji ') states respectively. For comparison, we note that the 
constraint arising in this case for the quantum bouncer is 

From the momentum-completeness sum rule in Eqn. fllSp . we use the evaluation of the 
the p"^ operator from Eqn. f lT5]) and find 



all A: 



VkiVk - Cn)^ 3 



while the corresponding sum rule constraint for the bouncer is 

For the sum rule in Eqn. ( 1201) . the fact that [V^'(-2)]^ = (±-^)^ = F^ for both positive and 
negative values of z, yields 
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15| . For that system, we 



The convergence of these sums is clear since the large k behavior of the Cfc ~ k"^^^ implies 
that each term scales as k~^'^. Thus, p = 2 is the smallest power which will lead to a 
convergent sum rule. For the case of the quantum bouncer, this sum rule did not converge. 
For the sum rule in Eqn. flT9|) we require the value of V"{z). Given the cusp in the 
definition of V{z), we find that V"{z) = 2F5{z) so that the right-hand side of Eqn. (TT9|) 
becomes 

^l^i'^Wr- (34) 

We can then use the results of Eqn. (J9]) and immediately see that 

and note that the form of the right-hand side of this sum-rule constraint varies in form 

between the odd and even states. The only other case of which we are aware where the 

forces times momentum sum rule gives a similar result, namely where the sum rule depends 

on the wavefunction at the origin, is for the Coulomb problem 

have V^V^(r) = — e(5(r) and the sum rule gives a non-zero result for 5- wave (/ = states) 

only. 

To evaluate the constraint arising from the 2;-completeness relation in Eqn. fl27|) . we 
require the expectation values of z^ in the even and odd states. These can be derived by 
using the integral in Eqn. ( 1A9I) to obtain 

{^lV\i^V) = / (^) and (^W|.Vi^^) = P' (^ + ^) (37) 

where we once again notice the difference in form between the results for the two parities. 
Using these results, we find the following sum rules 

T6(n) = y , \ ,, = ^ (38) 

^^^^^Ki)S(c^^?f"i- ^''^ 

The corresponding result for the quantum bouncer ^ has a non-zero 'diagonal' term, as in 
Eqn. ( l27l) . which leads to the result 

S^ = Y ,^ ^ ^^^ = ^ (40) 
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It has been stressed [2] that the form of the standard second-order perturbation theory- 
result, namely 

\{n\V{x)\k)\^ 



E^'=Y. 



^ fR(0)_^{0)^ 



(41) 



can also be thought of as an energy-weighted sum rule. Jackiw [16>] discussed energy- 
difference weighted sum rules, containing factors such as {E\^ — Ey^^ with negative values 
of p. This concept was applied in Ref. |J] to obtain another independent constraint on the 
S^in) for the quantum bouncer by considering the perturbing effect of an additional constant 
force (linear field) via the Stark effect. 

The Stark effect for the symmetric linear potential has recently been discussed 8| where 
a straightforward expansion of the exact eigenvalue condition was shown to led to closed- 
form expressions for the second-order energy shifts. If the perturbing potential is defined 
as Vi^z) = Fz, then the results for the Stark shifts for the odd- and even-parity states 
respectively are found to be 











-) 




Using 
theorj 


the 
r are 


expression 


in Eqn. (|4lj). 


we i 


ind 








Ei-'^^ = 


-4 


fF 








^(..) ^ 


-4 


fF 
\F 



9 \F " 



and E(+.2) = _r(^) Ei+K (42) 

we find the the second-order shifts given by perturbation 



^=^\ 2 



Sn 



-F^\ 2 



^ £a 



.allfc 



allfc 



1 



which gives two new constraints, 

Trin) : 

U.in) -- 



E 



VkiVk-CnY 



VniCk -VnY 



Kn 



(43) 

(44) 



allfc 



Tlk{r]k - Cn 



36 



^Jft 



(Cfc - VnV 



5rjn 
' 36 



(45) 
(46) 



which can be compared to the the quantum bouncer, 

1 



S,{n) = ^ 



k^n 



(Cfc - Cnf 



36 



(47) 



Our earlier discussion of sum rules for the quantum bouncer [4J] allowed for a systemat- 
ically calculable hierarchy of constraints on sums of inverse powers of {C,k — Cn)- This was 
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made possible by repeated use of the commutation relation [x'^.p] = iqhx''~^ to recursively 
obtain sums over matrix elements for higher and higher moments. For the symmetric poten- 
tial, where matrix elements of odd powers of z vanish, that connection is lost and we know 
of no strategy to generate all of the Tp{n) and Up{n) in a systematic way. We can, however, 
make use of the two relations 

'^{n\x''\k){k\x\n) = {n\x''+^\n) (48) 

all k 

Y^{Ek-En){n\x''\k){k\x\n) = q (^j {n\x''~^\n) (49) 

all k ^ ^ 

for q odd to generate an infinite number of new constraints on inverse powers of [rjn — Ck) 
similar to the ones derived already, requiring only the expectation values on the right-hand 
sides, which can in turn be evaluated using the recursion relations in Appendix [Bl 

The sum rules involving dipole matrix elements for the symmetric linear potential, using 
Eqns. iHM and ( IT5l) . give constraints on sums of inverse powers of {rik — Cn)- In order to 
obtain constraints on differences of just the zeros of Ai'{() separately, namely the (77^ — rjn), 
we need to consider matrix elements which are non-vanishing for even-even states. For 
example, because of the parity constraints in the symmetric linear potential, the so-called 
monopole sum rule [l7| 

TiEk - E^)\{n\z'\k)\' = —{n\z'\n) (50) 

^ m 

kjtn 

will connect only odd-odd and even-even states; in this case the constraint that k ^ n is 
indeed required. For the odd-states, the corresponding constraint on the (n was found in 
Ref. y to be 



For the even-even case, we can use the new matrix element result in Eqn. (IA13P to evaluate 



{^i-v\^':') = - i!i^;+^^\^ (52) 



y- {Vn + Vkf ^ ]_ fSr]l J_ ^ 



which gives a more complicated constraint, namely 

k-tn '^n^kiVn - rjkY 36 V 15 ' 5r]n 
Once again, we have checked all of these results numerically using Mathematical 
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The new classes of constraints on the zeros of Ai and Ai' derived here, from apphcation 
of sum rules for the symmetric linear potential, are seen to be qualitatively similar to those 
obtained from the parity-restricted version of this potential, namely the quantum bouncer. 
In contrast, the nature of the mathematical relations dictated by sum rules applied to the 
parity- restricted version of the harmonic oscillator are qualitatively very different, as we will 
see in Sec. |Vl after first reviewing sum rules in the standard oscillator system. 

IV. REVIEW OF SUM RULES FOR THE HARMONIC OSCILLATOR 

Before examining the parity-restricted version of the harmonic oscillator, we briefly review 
the solutions for the the familiar oscillator potential, as well as the structure of the quantum 
mechanical sum rules for this system. The solutions for the Schrodinger equation 



2m dx^ 2 

can be written in the form 



^„(x) = -^Hn{y) e ^ /^ where c„ = , "" , ^ (55) 

V/3 V2"r 



"■nlwir 



where 



h 

X = /3y and P = \l (56) 

V mco 

with y dimensionless and the Hn{y) are the Hermite polynomials. The corresponding energy 
eigenvalues are 

2E 
En = {n + l/2)huj or e„ = — ^ = 2n + 1 (57) 

with n = 0, 1, ... so that the differential equation in dimensionless form can be written as 

C(y) = (y'-en)V^n. (58) 

The solutions have parity given by P„ = (—1)" and the expectation values of any odd power 
of X vanishes, so that, for example, (n|x^^"'"^|n) = 0. 

Both X and p can be written in terms of raising and lowering operators, A and A\ and 
hence the multipole matrix elements exhibit an exceptionally simple structure leading to 
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absolute selection rules. For example, we have from standard textbooks the relations 

{n\x\k) = —= l6n,k-i'/k + 6n,k+Wk + l\ (59) 

(n|x2|fc) = ^ ^5n,k~2^k{k-l) + 5n,k{'^k + 1) + ^n,k + 2^{k + l){k + 2)^ (60) 



{n\x''\k) = ^ {5n,k-WKk -l){k-2)+ 35„,fc_ifc3/2 + 5„ ,^,(A; + 1)3/2 

+5n,fc+3V(^ + l)(A; + 2)(A; + 3)} . (61) 

These, and similar expressions for matrix elements of powers p, imply that all of the familiar 
sum rules discussed above will be 'super-convergent', namely that the relevant infinite sums 
will actually be saturated by a finite number of terms, and hence satisfied in a trivial way. 
We will find that the situation is dramatically different in the case of the oscillator restricted 
to the half-line, at least for the matrix elements of odd values of x, and that is the subject 
of the next section. The matrix element relations for even powers of x, including Eqn. ( 160|) 
will still be relevant for the restricted oscillator case with minor relabeling. 



V. THE PARITY-RESTRICTED HARMONIC OSCILLATOR 

The solutions to the quantum mechanical problem of a particle in the potential in Eqn. ([3]) 

are easily obtained from the odd-parity solutions in Eqn. f l55|l . In dimensionless notation, 

we have 

r,, /v^^2n+i(l/) fori/>0 

^n[y) = < (62) 

[^ for y < 

where n = 0,1,2,... for the 'half'-SHO states ipn{y) associated with the odd solutions of 
the oscillator, and an appropriate change in overall normalization. We will henceforth use 
integrals over either the dimensional x or dimensionless y variable as deemed most useful. 
We note for future reference that the derivatives of the solutions at the origin necessary for 
matrix element calculations are given by 



m - V2c,„« Hi„,M ^ ^J;;';;^ (2n + DH ^ ^J^^^ (63) 



and 



[^-(0)f= ' (2n + l)!!^^ 
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■4(2n + l) 
22"fn!)2 



t (^^) 



where 

4(2ra + l)! 

22"fn 



A.-^SS^- m 



Using the Stirhng approximation, n\ ~ y/27m{n/e)"', we have for large n 



D,,^8J- (66) 

which will be useful in examining the semi-classical limit. 

The quantized energies for the 'half'-SHO are then given by En = i?2n+i = (2ri + 3/2)hu 
or e„ = 4n + 3. The ipniv) form an orthogonal set which can be seen explicitly by using 
the recursion relation derived for oscillator solutions in Appendix [Bl Using f{y) = 1 in 
Eqn. (JBS]), we find that 

{en-emy / -ipniy) V^m(y) dy = (67) 



^0 

so that (ipmli'n) = if n 7^ m. 

For the important dipole matrix elements, we use f{y) = y, so that f'{0) = 1 and find 
that 

{^m\ym = ryV^n(y)V^n.(y)^l/= - ^.5^°^"^":^ ^ .1 (68) 

Jo 2[4(n - my - 1\ 

since e„ = 4n + 3. 

For the special case oi n = m, the expectation values in the state ipn are therefore given 
by 

(V^nlxl^.) = -^/3 -^ ^/3 (69) 

in the large n limit. We can confirm this using the classical (WKB-like) probability distri- 
bution for the 'half'-SHO, namely 

P^?l{x) = -^=^= where < x < A^ (70) 

where the upper classical turning point, A^,, is given by 

-muj'^Al = En = nw{2n + 3/2) or A„ = /3V4n + 3 . (71) 

The classical expectation value is then 

inMnW^ rxPl;i(x)i.^?^ = ?^^l±Il ^ fV«V (72) 

Jo IT n V ^ / 



for large n, which agrees with Eqn. fl69|) . 
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The dipole matrix elements in Eqn. (!68|) are clearly very different from the 'full'-SHO case, 
with each state being connected to all of the others in a very simple form, proportional to the 
derivative of the wavefunction at the origin for each state, and with an 'energy denominator' 
factor. This is the only example of a dipole matrix element in a model ID system of which 
we are aware for which the diagonal (n = m) and off-diagonal {n ^ m) cases can be written 
with the same simple expression. The corresponding sum rules will then be realized in a 
completely different manner than the super- convergent form for the more familiar oscillator. 

For example, the TRK sum rule in Eqn. (IT71) gives the constraint 

>^ {k-n)Dk . . 



where the Dn are given in Eqn. ( 1651) . In a similar way, the x-completeness relation of 
Eqn. (EI]) gives 

where we use the fact that 

{i>n\y^\i>n) = i2n + 3/2) (75) 

making use of the matrix element in Eqn. ( 160|) with a suitable relabeling. One can then 
combine Eqns. f l73|) and f l75|) to obtain the constraint 

^"E r4(^_^MLii2 = (4^ + l)(2^+l)^- (76) 

Looking at the n and k dependence of the dipole matrix elements, one can confirm that the 
sum rules in Eqns. (TT9|) and (!20|) which depend on derivatives of V{x) are not convergent, 
this result is consistent with the nature of the discontinuous potential. 

Matrix elements involving even powers of x can be written in terms of the results for the 
'full'-SHO and have the same simple 'nearby neighbor' selection rule structure. For example, 
a simple relabeling of Eqn. fl60|) gives 



i^nly^l^k) = \ [5n^k-W{2k + l){2k) + 5„,fc(4A; + 3) + 5^^u+W {2k + 2){2k + i)'^ (77) 

and the monopole sum rule in Eqn. f l50|) is saturated by a finite number of terms. Ma- 
trix element relations using odd powers of x, such as those in Eqns. ( HHj) and (H9|) give 
increasing complicated relations involving the D^ n since the matrix elements of x^'^"'"^ all are 
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proportional to {il)n\x\il)k) ■ For example, using the recursion relation in Eqn. (IBSp . we find 
that 

^ -3(2n + 2A: + 3) (^<(0)^\ 

[4(n - A;)2 - 9] [4(n - A;)2 - 1] (^ 2 j ^ ' 

One can continue to generate increasingly complex constraints by use of Eqns. ( I48l) and (H9l) 
for the non-trivial case of matrix elements of odd powers of x. 

The realization of the sum-rules for the parity-restricted version of the oscillator are 
completely different than the trivial manner in which they are satisfied for the ordinary 
oscillator system and generate an infinite number of constraints on the Dn- 

VI. CLASSICAL VERSUS QUANTUM MECHANICAL RESULTS 

We have seen that some of the matrix element expressions and/or sum rules for the 
symmetric linear potential yield identical results when considering both the even and odd 
cases, with only the substitution ^„ <-)■ ?]„ required: the virial theorem result in Eqn. ( 1T2ll and 
the TRK, momentum completeness, and force-squared sum rules in Eqns. (128!) . (!30l) . (1311) . 
and fl33|) . For others, however, such as the second-moment expectation values in Eqn. ( !37|) 
or the force times momentum result one finds slightly different results. 

To examine these small differences, we will consider the classical underpinnings of some of 
the results which arise in the evaluation of expectation values. By generalizing the recursion 
relation derived by Goodmanson [10|] in Sec. [B]one can generate the expectation values and 
off-diagonal matrix elements of any power of z. 

For example, from Ref. (^ we know that expectation values of z'^ for the quantum bouncer 
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solutions are given by 

{n\y\n) = ^ (79) 

{n\y^\n) = ^ (80) 

{n\y'\n) = ^ + ^ (81) 

{n\y'\n) = ^^ + ^ (82) 

^ '^ ' ^ 315 63 ^ ^ 

, , 5, , 25QC 1808C ,oo^ 

(n y^ n) = — H ^ (83) 

^ '^ ' ^ 693 3003 ^ ^ 

where y = z/p and we only show the dimensionless results. Each result has a highest-order 

term of order (Cn)^, followed by sub-leading terms of order Cn~^'^5 if present at all. Using the 

large-n expansion of the (n in Eqn. ([7]), we see that the sub- leading terms are a factor of 

{Cn)~^'' ~ n~'' smaller and so become negligible in the classical limit. This suggests that the 

leading terms are indeed what one would expect from a purely classical probability density. 

To confirm this, we note that for the quantum bouncer we have 

P^"2 (^) = , ^ (84) 

where An is the upper classical turning point, defined by E^ = FA^. If we equate the total 
energy with the quantum mechanical result En = {Fp)(„ and write An = p(n, the classical 
probability density reduces to 

P^"2 (z) = , ^ . (85) 

2pVCn(Cn - z/p) 

We briefly discuss, in Sec. [Cl how this classical distribution can also be extracted directly 
from the large n limit of the exact quantum solutions. 

The expectation values of moments of position are then given by 

{n\zP\n)cL = / zPP^^l{z)dz 
Jo 

= pP / — dz 

Jo 2^yCniCn - Z) 



(86) 



2 Jo v^T^ 
^5(P+ 1,1/2) 

(pCn)^r(i+p)r(i/2) 

2 r(p + 3/2) 
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and this expression agrees with the leading order terms in Eqns. ( !79|) - ( 183|) up through p = 5. 
Using the recursion relation of Goodmanson, reviewed and extended in the Appendix, we 
can assume generally that the behavior of the leading term for the expectation values can 
be written as {n\x^\n) = A^ and the recursion relation requires that 

which agrees with the classical result in Eqn. fISB]) for all p values. 

For the symmetric linear potential, the expectation values for odd powers of x vanish, 
but the integrals in Eqns. f lTU]) - dHSD are still useful. For example, for the expectation value 
of V{z) = F\z\, we require the result in Eqn. fl79|) . Moreover, because of the piecewise 
continuous definition of the wavefunctions in Eqns. dl]) and ([6]), the integrals obtained by 
using the recursion relations in Appendix |B] are necessarily defined over the interval (0, oo) 
and then extended over all space, so intermediate results for integrals over the half-line which 
eventually vanish due to parity constraints can still be useful. 

To compare the leading and sub-leading contributions to integrals used for the expectation 
values for odd and even states, we can compare the results in Eqns. (!79|) - (!83|) to similar 
results for the even states. These can be defined by using the states in Eqn. ([6]), integrated 
over positive values of x and normalized as for the quantum bouncer. For those cases we 
find 

{n\y\n) = ^ (88) 

H.^|.)^f + ^ (89) 

{n\y'\n) = ^-^ + ^- (90) 

, , 4, , 128r7^ 64r]n ,„,, 

(n\y^\n) = - + — - (91) 

^ '^ ' ^ 315 45 ^ ^ 

(nly^ln) = H H (92) 

We note that the leading terms in each case (and quite generally for all values of p, using 
a recursion relation argument as above) are identical (with C„ -H- ?7„), but that the next-to- 
leading orders reflect differences between the classical and quantum probability densities. 
These differences vanish in the large n limit. 

In the context of the 'half '-SHO, we have already noted that there can be similar agree- 
ment between the exact quantum-mechanical expectation value of x in Eqn. fl69|) and the 
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corresponding classical result in Eqn. ( !72|) in the large n limit. One can use the recursion 
relations in Eqn. (IB 8 1) to obtain the highest-order terms in the expectation values of x^ 
and one again finds agreement with the semi-classical results for large n. For example, the 
classical result is 

/ I PI \ 2 [^- x^ ^ 2{A^Y f' jfdy {A^y T{{l+p)/2) 

{n\x' \n)cL = — I — dx = / — = — ^ — -— —— . 93 

Given that traditional sum rules involve transition matrix elements, there is no reason to 
expect that semi-classical probability arguments will provide any useful information on their 
evaluation. One important exception, however, is the form of the second-order perturbation 
theory result for the energy, as in Eqn. (HTj) . which is of the form of an energy-difference 
weighted sum rule. In that special case, semi-classical expressions for the quantized en- 
ergy, such as the WKB approximation, can sometimes provide guidance on the form of the 
energies, at least in the large n limit. 

An example of such a connection is the use of approximate WKB-type methods in the 
evaluation of first-order perturbation theory results using classical probability densities, as 



in Ref. 18|. More surprisingly, it has been pointed out that WKB energy quantization 



methods can give the correct large n behavior of the second-order energy shift due to the 



Stark effect in two familiar model systems, the harmonic oscillator and infinite well 19|. 
This approach was used in Ref. (Sj where the exact result for the second-order energy shift 
due to a constant external field for the symmetric linear potential was derived for the first 
time, giving the results in Eqn. fH2|) . The WKB prediction for this case is given by the 
quantization condition 

V2^ \J En - {F\z\ + Tz) dz = {n + l/2)hTT (94) 

where n = 0, 1,2, .. and the classical turning points are A± = ±En/{F ± F). The WKB 
prediction for the energies are then 

where En = £o{3n(n + 3/4) /2)^/'^ is the zero- field WKB approximation, which agrees with 
the exact results in Eqn. ([7]) for large n. This implies that the first-order Stark shift vanishes 
(as it must by symmetry) and that the second-order terms are 

4=' = -5(?)'^"°'- (^^' 
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As pointed out in Ref. J8|, this semi-classical result brackets the two exact quantum me- 
chanical expressions for the even and odd states in Eqn. P2|) . giving it as the 'average' 
effect. 

Prompted by this partial success, we wish to examine to what extent a similar WKB-type 
analysis will give reliable answers for the first- and second-order energy shifts due to the 
Stark effect for the 'half '-SHO discussed above. We first note that the WKB result for the 
energy eigenvalues for the potential in Eqn. ([3]) without an external field are given by 

y/^m / ^/En - muj^x'^/2 dx = {n + CL + CR)h7r (97) 

Jo 

where Cl, Cr are the appropriate matching constants for an infinite wall and 'linear' poten- 
tial respectively, given by Cl = 1/2 and Cr = 1/4, respectively. Evaluating the integral, 
the WKB prediction is then given by En = {2n + 3/2)hw, reproducing the exact result. 
The corresponding expression including a perturbing linear field, V{x) = Fx, is then 

^ f^""' I — 

\r2im \ \IE - muj'^x'^/2 -Fxdx={n + ?,/A)n-K (98) 

Jo 
where the upper turning point is given by energy conservation to be 



y2mu^En + F^ - F 
Al^^ = ^ -, • (99) 



mu'^ 



The integral can be done in closed form, but we only require the result expanded to second 
order in F, which gives 



Tc / m 



2 V Tau'^ 



2 I2E„^ F^ 



En--\ ^^F + 



{n + 3/4)h7i. (100) 



TT V mu'^ 2muj'^ 

This can be rationalized to give the simple quadratic equation 

En = Rn± VK. - Zl (101) 

where 

f' 4f' 
Zn = - ^ and Rn = Zn + - -. (102) 

Solving for En-, again as a series in F, we find that 

E^^\WKB) = {2n + 3/2)huj (103) 

E(^\WKB) = ^./^ (104) 

TT V "fT^^ 

-^ + 4) A- (105) 

2 IT'' J mu^ 
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The zero-order result is the standard WKB prediction noted above, while the first-order ex- 
pression coincides with first-order perturbation theory (using the matrix element in Eqn. f l69|) 
in the large n limit. 

The second-order result is more interesting. The exact quantum mechanical result for the 
second-order Stark effect for the ordinary harmonic oscillator is En = —F /2muj'^ which is 
most easily obtained by a simple change of variables in the original Schrodinger equation, 
and trivially confirmed in second-order perturbation theory, using the matrix elements in 
Eqn. (l59l) . For the 'half'-SHO, the WKB prediction is still a constant negative shift, the 
same for all states, but with a non-trivially different coefficient. We can find no simple way to 
extract the exact second-order result from a direct solution of the differential equation, but 
the second-order perturbation theory expression in Eqn. ( 14T1) can be evaluated numerically 
using the dipole matrix elements in Eqn. (1^5]) . 

To compare the first- and second-order predictions from the WKB approach with the 
exact results from first- and second-order perturbation theory (PT), we plot the differences 
between the WKB and PT methods in Fig. [H as a function of the quantum number n. As 
mentioned above, the first order predictions in both approaches agree in the large n limit, 
which we've demonstrated here analytically. We note the more interesting result that the 
simple expression in Eqn. f llOSp gives the correct large n behavior of the second-order Stark 
shifts for the 'half'-SHO problem, reproducing a non-trivial numerical factor which would 
otherwise be difficult to extract. 

We note that this is another example of where a WKB approach to the evaluation of 
the second-order Stark effect correctly predicts the large n behavior. This further justifies 






the discussion Ref. llSl where such an approach was used to systematically evaluate the 
second-order energy shifts due to an external field for general power- law potentials, Vk{x) = 

Vo\x/a\''. 

VII. CONCLUSIONS AND DISCUSSION 

We have examined two cases of parity-related potentials, the quantum bouncer extended 
to the symmetric linear potential, and the harmonic oscillator reduced to the 'half'-SHO, 
in order to probe the importance of the continuity of the potential on the convergence 
of quantum mechanical sum rules. We have indeed seen that the smoothness of V{x) has 
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clear consequences on which sum rules will be realized, which in turn is closely related to the 
convergence properties of the individual matrix elements. For the symmetric linear potential, 
we find new constraints on the zeros of the derivative of the Airy function, but note that 
they are very similar in functional form to those derived from the quantum bouncer. On the 
other hand, the infinite set of constraints which arise from the 'half '-SHO are qualitatively 
very different from the super-convergent sums found in the realization of sum rules for the 
more familiar oscillator. 

The study of parity-related potentials is also motivated by the desire to find examples 
where there is a substantial overlap between the energies and wavefunctions (needed in 
the evaluation of matrix elements) of two quantum systems as they relate to sum rules. 
That connection is realized here by the fact that the odd states of a symmetric potential 
remain solutions of the parity-restricted version, so that all of the resulting energies and 
wavefunctions (save a trivial normalization) are still solutions of the parity-restricted partner 
potential. 

Another class of quantum mechanical problems which has similar strong connections 
between the energy levels and wavefunctions are super-partner potentials 20| , V^'-"-' (x) and 
V^~^\x), in the context of supersymmetric quantum mechanics (SUSY-QM). In that case, the 
spectra of the two systems are identical, except for the zero-energy ground state (Eq = 0) 
of the first system, which is absent in the second. Another reason for us to consider the 
'half '-SHO potential in such detail is that it is easily extended to generate an appropriate 
V^'^\x) pair in SUSY-QM. For example, the potential 

1 3 

V^-\x) = -mu'^x^ - -hw fora;>0 (106) 

has the energy spectrum En = 2nhuj for n = 0, 1, 2, ... with the ground state wavefunction 

Mx) = /^ e-("/^)'/2 for X > (107) 



and the remaining states given by Eqn (1621) . 

Using this ground-state solution to form the super-potential, we find 

y/2^\Mx)J V2^\x (3^ J ' 

allowing us to construct the super-partner potential 

1 iW' 1 

\/(+)(x) = -mu;V + ---hjj. (109) 
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This has the form of the radial equation for the three-dimensional harmonic oscillator with 
the special choice of the angular momentum quantum number / = 1 (giving the /(/ + 1) = 2 
factor in the centrifugal barrier term) and an overall constant shift in energy. Using standard 
results for the energy eigenvalues for that system, we find that 



EW =fi^(2n + l + -j --hw = 2nw{n + l) 



(110) 



for the relevant / = 1 case. One can also use standard textbook results to obtain the properly 
normalized solutions to be 



i,^^\x) = iV„x^e-(^/«^/^L(3/2)(xV/3^) with iV„ = J ^5(2fc + 3)n^ - ^^^^^ 



This is one example of many familiar super-partner potentials 20| which can be systemati- 
cally studied in the context of quantum mechanics to probe the delicate interplay between 
energy level differences and matrix elements which must exist to guarantee the realization 
of the infinite number of sum rules which one can generate using the simple procedures 
outlined here. 
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Appendix A: Necessary Integrals Involving Airy function 

In this Appendix, we collect indefinite integrals of products of Airy functions, useful for 
normalizations, expectation values, and matrix elements as first noted by Gordon [l3| and 
Albright [14], adding one new result (Eqn. ( IA13I) . 



We then extend the recursion relation for general matrix elements of Airy functions, first 
derived by Goodmanson [10| , necessary for both the even and odd states of the symmetric 
linear potential in Sec. [B] and also generalize those recursion relations to the case of the 
harmonic oscillator problem, allowing one to evaluate matrix elements for the 'half'-SHO 
problem. 
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We first assume two arbitrary linear combinations of the Airy difi^erential equation, 

A{C -(3) = aAziC -/3) + bBz{C - (3) and S(C - P) = cAi{C - P) + dBz{C - P) (Al) 
where 

^"(C - /3) = (C - miC - /3) and B"{C -/?) = ((- P)B{C - /?) • (A2) 

In the expressions below, when we refer to integrals of such solutions corresponding to 
different values of the 'shift', /3i 7^ /32, we write Ai for A{C — /3i) and B-^ for B{C — (32)- 

Each of the indefinite integrals we need can be written in terms of four basic quantities, 
namely 

Fi = {AB} (A3) 

F2 = {A'B-AB'} (A4) 

F3 = {A'B + AB'} (A5) 

F4 = {A'B'} (A6) 

times polynomials in ( and the appropriate values of /3. 

For the integrals involving solutions with the same 'shift', we have 



ABdC = (C-/3)Fi-F4 (A7) 

I (ABdC = i(c2 + /3C-2/32)Fi + lF3-i(C + 2/3)F4 (A8) 

CABdC = 4(3C' + /3C' + 4/3'C- 8/3=^-3)^1 + 4(30 + 2/3)^3 
15 15 

-4(3C' + 4/3C + 8/32)F4. (A9) 
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For integrals involving derivatives, we can modify one of the identities in Albright IJ] to 
write for shifted solutions 

JA'B'dC = _i(C-/3)2Fi + ^F3 + i(C-/3)F4. (AlO) 

Finally, for integrals with two different shifts (/3i 7^ P2), "we have results for p = 0,1 
from Gordon [13[ and a new result derived here for p = 2. Recall that the T, below will be 
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constructed from the terms in Eqns. (]A3P - (]A6|) using Ai = A{( — f3i) and B2 = B{( — (32). 

1 



Ai B2 dC 



(/32-/3i) 



(All) 






C 2 

+ 



(/32-/3l) (/32-A)3 



CA,B2dC 



(A - /32)^ 

12(/3i + /32) , 2(-12+(/3i + /32)(/3i-/32)') 



+ 



L(/3i-/32)^ (/3i-/32)^ 

4(-6 + (/3i + /3i)(/3i-/32)') 12 



C 



(/3i-/32)^ 



+ 



(/3i - /32)^ 



(/3i - P2f 



F.+ 



rC- 



1 
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L(/3i-/32)^ (/3i-/32) 



(/3i-/32)3' (/3i-/32) 



C^ 



+ 



(A12) 

Fi 

F2 
(A13) 



Appendix B: Recursion relation for matrix elements of Airy functions and harmonic 
oscillator solutions 



Goodmanson 



10| derived a recursion relation relating the power-law matrix elements of 



solutions for the quantum bouncer problem, involving integrals of the form 

/■oo 

R,= / CAi{C-Cn)Ai{C-Cm)dC (Bl) 

Jo 

assuming that the solutions satisfied the boundary condition for the bouncer problem, 
namely that Ai{—Cn,m) = 0. One can easily repeat his analysis, not making that assump- 
tion, to find a more general recursion relation valid for integrals relevant for both the even 
and odd states of the symmetric linear potential. We can also repeat the analysis to obtain 
relations on matrix elements of the solutions of the harmonic oscillator, which will prove 
useful for the parity-restricted 'half '-SHO in Sec. El 

a. Recursion relations for Airy functions 



We begin by adopting the notation 

An = AniC) = Ai{( - (n) whcrc 



/l_ 



<(C) = (C - Cn)A. 



(B2) 



and then assume a general f{() as a well-behaved function, at least multiply differentiable 
and satisfying lim^^oo /(C)^^[C]^ = 0. We will, in fact, use f{() = ('^ to evaluate various 
matrix elements. 
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We start with the mathematical identity 

and then use integration by parts muhiple times to isolate the An{C)Am{C) terms times 
derivatives of f{() (on the left hand side) and resulting surface terms (on the right hand 
side.). The resulting integral can be written in the form 

An An. [/(»)(C) - 4(C - Ca.e)r(C) " V\0 + (Cn - Cmf f [O] ^C (B4) 

which can be shown to be equal to 

/ = WiO^^A'^ - nO{AnAm)']7 + [nOArAm]'^ ' [2A„A„(C " Ca.e)/'(C)]o~ 

-[(c„-u«^™-^nAj/(c)r 
= -2/'(ox(o)a:„(o) + r (0) K(o)A„(o) + Ar^ALm 

-r'(0)A„(0)A„(0) - 2A„(0)A„(0)f (0)Ca.e 

+ (Cn - Cn^)/(0) K(0)A„(0) - A^A^m (B5) 

where C^ave = (Cn + Cm)/2. The first term on the right-hand-side reproduces the Goodmanson 



10| result. We have assumed that f{() is well-enough behaved that surface terms at infinity 
vanish. One can then evaluate matrix elements over the range (0, oo) recursively for any 
combination of products of either even or odd solutions, where simplifications occur since 
A'nmi^) = or An,m{0) = respectively. 

b. Recursion relations for oscillator solutions 

For the ordinary harmonic oscillator, expectation values and matrix elements will be 
evaluated by integrals of solutions over the range (—00, +00). But an almost identical 
derivation to that above allows for the evaluation of related quantities over the half-line 
(0, +00) which will be important for the parity-restricted oscillator (the 'half'-SHO). For 
example, if we label ipn{y) as a solution to the dimensionless oscillator problem in Eqn. (l58ll . 
namely 

€(2/) = (2/'-en)^n, (B6) 

we can make repeated use of that relation and integrations by part to write 

/■oo 

J ^ / My) ij„,{y) [/(-)(2/) - 4(y2 _ ea.e)f"{y) - Ayf'{y) + (e„ - em)'f{y)] dy (B7) 
Jo 
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in terms of surface terms. In this way, we find that 

J = -2/'(0)^;(0)V^:„(0) + V^„(0)V'^(0)[r(o) + 2/'(0)e,,e] 

r(o)K(o)^^(o)+^„(o)^:„(o)] 

-(en - e^)f{0) K(0)^™(0) - MO)^'mm (B8) 

where eave = (cn + em)/2. Note that for the 'half '-SHO, only the first term where '?/'^(0) and 
ipmi^) are non-zero 

We can confirm that this expression encodes the absolute selection rules for matrix ele- 
ments of the standard oscillator, as in Eqns. (13^ - (IMj) . For example, for the dipole matrix 
elements, we know that for the oscillator only even-odd matrix elements will be relevant, so 
using f{y) = y, we find that 

J = [{en - eraY - A] ^^iv) ^miv) dy = f) (B9) 

so that the dipole matrix elements vanish unless {n — rnf' = 1 oi n = m ± 1. 

Appendix C: Semi-classical probability density for the quantum bouncer 

It is straightforward to see how the classical probability distribution in Eqn. (18^ for 
the quantum bouncer arises as the locally averaged limit of the exact quantum mechanical 
solution. For the bouncer, the solution is given by 

^„(z) = V2iv(-)A^f--cl (CI) 



where Nn is given by Eqn. (jHI). Standard handbook HI] results exist for the behavior of 
Ai{C) and Ai'{C,) for large negative values of C, namely 

A^{-0 - ^-n-^/^in (e + 9 (C2) 

A^{-0 ~7r-VV^/^cos(e + ^) (C3) 

where ^ = 2a;^'^/3. The quantum- mechanical probability density is then given by 



P^^liz) = \Uz)\' = - 



Ai{z/p-Cn) 
A^'{-Cn) 



(C4) 



^ sin^ [2(Cn - z/pfl^/^ + vr/4] ^^^^ 



P\/Cn{.Cn - z/p) 
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where we use use the expression for (n in the large n hmit in Eqn. ([7]) to write 

2/3\ 3/2 



cos(^e + ^)~cos^ 



Y(--1/4) 



+ j}= cos{nn) = ±1 (C6) 



In the semi-classical limit, the oscillatory component locally averages to 1/2, giving the 
purely classical result in Eqn. ( l85l) . This explicit derivation agrees with the standard result 
for a purely classical probability distribution for position, namely 



^-(^) = 7\/w^TM) (c^) 



where r is the classical period, given by 

"'' dx 



(C8) 



2 Ja Vi^) 

with v{x) being the classical speed and a,b are the classical turning points for bounded 
motion. 
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FIG. 1: Comparison of the WKB estimate (WKB) of the first- and second-order energies for the 
Stark effect (from Eqns. (jl04p and (|105p ) versus the exact results from perturbation theory (PT). 

(1 2) il 2) 

The quantity En ' (WKB) /En ' (PT) — 1 is shown as dashed (first-order) and dashed (second- 
order) respectively. 
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